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Light diffusion is usually associated with thick, opaque media. Indeed, multiple scattering is 
necessary for the onset of the diffusive regime and such condition is generally not met in almost 
transparent media. Nonetheless, at long enough times, transport in an infinite thin slab is still 
determined by a multiple scattering process whose complete characterization is lacking. In this 
paper we show that, after a short transient, the mean square width of the transmitted intensity 
still exhibits a simple linear increase with time as predicted by diffusion theory, even at optical 
thickness as low as one. Interestingly, such linear growth is predicted not to depend neither on the 
slab thickness nor on its refractive index contrast, yet the accuracy of this simple approximation in 
the ballistic-to-diffusive regime hasn’t been investigated so far. By means of Monte Carlo simulations, 
we find clear evidence that boundary conditions play an active role in redefining the very asymptotic 
value of the diffusion coefficient by directly modifying the statistical distributions underlying light 
transport in such media. In this respect, we demonstrate the need to distinguish between a set of 
intrinsic and effective transport parameters, whose relation and interplay with boundary conditions 
remains to be fully understood. 

PACS numbers: 42.25.Dd, 87.64.Aa, 87.64.Cc, 42.62.Be 


I. INTRODUCTION 

The modeling of light transport in disordered media 
represents an intensive research field both for its funda¬ 
mental interest and for its applications. In this ubiqui¬ 
tous class of media, because of the stochastic nature of 
multiple scattering and absorption, light transport can 
effectively be described as a random walk of energy pack¬ 
ets mm- Whenever an unbounded or thick geometry 
can be assumed, the diffusion approximation (DA) pro¬ 
vides an incredibly robust and useful simplification of 
such problem since it casts a complete set of analytic 
predictions on all transport observables, also in the slab 
geometry [3]. Their validity relies on the assumption 
of almost isotropic radiance, which is underlying in the 
diffusion approximation but progressively breaks down if 
the characteristic scattering length becomes comparable 
to absorption or the thickness of the slab. In the latter 
case a significant fraction of light is usually able to go 
through the sample without undergoing any scattering 
event and the slab becomes semitransparent. This con¬ 
figuration is particularly relevant in many applications, 
ranging from coatings to photovoltaic and biomedical, 
and typically requires a refined analysis to overcome the 
limitations of the diffusive approximation. In this case 
Monte Carlo (MC) simulations are typically considered 
as the gold standard method to solve the radiative trans¬ 
fer problem, in that they allow a direct implementation 
of the random walk model with accurate boundary and 
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geometric conditions [4]. With the continuous growth of 
computing power and parallelization capabilities, Monte 
Carlo simulations are gaining momentum as a convenient 
solver both for the forward and inverse light transport 
problem. Indeed, by resorting to MC simulations, novel 
insights are continuously unveiled on the peculiar physics 
that governs light transport in complex media such as het¬ 
erogeneous or anisotropic materials 011]. Conversely, we 
show that our comprehension of way simpler systems such 
a homogeneous, isotropic single slab is still incomplete. 

To this purpose we investigate by means of Monte Carlo 
simulations what is usually referred to as the diffusive- 
to-ballistic transition in the thin slab geometry, with a 
particular focus on transverse transport. Indeed the thin 
slab geometry is typically associated with a ballistic trans¬ 
port regime and consequently to the inapplicability of 
diffusion theory. Instead, the former does not necessarily 
imply the latter, which becomes clear when addressing 
in-plane transport where a deeply multiple scattering 
regime can indeed occur [7j. In an infinitely extended 
slab, transverse propagation is unbounded and will even¬ 
tually become diffusive at sufficiently long times. Nonethe¬ 
less, even if asymptotic, we will show that this diffusive 
regime is not the one predicted by diffusion theory based 
on the intrinsic transport properties of the investigated 
medium. In this respect our findings notably show that, 
for example, even a minute tuning in the refractive index 
contrast allows to sensibly modify the long time behavior 
of semitransparent slabs, which might therefore be still 
functionalized in terms of scattering properties in spite of 
their low optical thickness. 

The paper is organized as follows: in Section |ll] it is 
shown that a typical transverse diffusion signature such 
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as the linear growth of the mean square width (MSW) 
is still clearly visible after a short transient even at an 
optical thickness as low as 1. In Section m we find and 
characterize deviations of both the expected mean square 
width slope and transmittance decay lifetime from their 
expected values according to the diffusive approximation, 
discussing also experimental implications. In Section [IV| 
we describe the origin of the observed deviations in terms 
of a statistical analysis of simulated transport, unveiling 
the peculiar features of a diffusion process which remark¬ 
ably differs from the one expected by standard diffusion 
theory. 


II. TRANSVERSE TRANSPORT IN A SLAB 
GEOMETRY 


The slab geometry (see Fig. is a very simple yet 
relevant model in a wide range of applications. As such, 
transport properties in this configuration have been ex¬ 
tensively addressed in the literature and tested against 
diffusion theory, but mainly focusing on axial rather than 
lateral migration [EHH]. This is surprising since the the¬ 
ory fully describes transverse transport as well, providing 
the simplest relations linking experimental observables 
with the transport properties of the sample. According to 
the standard diffusion approximation, the instantaneous 
spatial profile emerging from a slab characterized by a 
diffusion coefficient D is always expected to be a Gaussian 
distribution I{p, t) (x exp(—with a mean square 
width growing as w'^ {t) = 4Dt. Such width represents the 
variance, which is more generally defined for an arbitrary 
distribution I (p, t) as 


w^it) 


j p^I{p,t) dp 
//(p,t)dp 


( 1 ) 


The mean square width value is most remarkably inde¬ 
pendent of absorption (which cancels out exactly at any 
time m) and, according to the simple diffusive model, 
independent of both the slab thickness and its refractive 
index contrast with the surrounding environment. This 
stands in great contrast with any other expression avail¬ 
able for the slab geometry in the diffusive approximation, 
where the correct assessment of absorption, thickness and 
extrapolated boundary conditions play a critical role in 
correctly retrieving the optical properties of the sample 
nmn], thus representing a largely appealing quantity 
to be experimentally investigated [TH]. Notably, even at 
an optical thickness as low as 1, a perfectly linear mean 
square width growth is permanently recovered soon after 
a short superlinear transient, as Figureillustrates for an 
index matched sample with unitary refractive index. As 
shown in the following, a qualitatively similar behavior is 
actually common to a broad range of optical parameters, 
which further makes the mean square width a valuable 
experimental observable also at low optical thicknesses. 


III. MEAN SQUARE WIDTH EXPANSION 


The validity range of the simple w'^{t) = ADt predic¬ 
tion has never been tested, to our knowledge. A compre¬ 
hensive Monte Carlo investigation, whose complete set 
of results are reported in a related work m, revealed 
the presence of non-trivial discrepancies which challenge 
the established framework linking radiative transfer to 
diffusion theory. For a single slab geometry with fixed 
thickness Lq = 1000 pm we extracted the transmittance 
mean square width growth rate from a large set of dif¬ 
ferent simulations in the (n, p, OT~^) parameter space, 
where n = rijn/riout is the refractive index contrast, g is 
the scattering anisotropy factor and OT“^ = k/Lo — a 
dimensionless measure for the transport mean free path 
It — is the inverse of the optical thickness. For each 
simulation we estimate D by performing a linear fit on 
the mean square width w‘^(t > At) as a function of time, 
where t is the simulated transmittance decay lifetime 
and At is selected as a lower limit so to always exclude 
early-time energy packets transmitted before the onset of 
the diffusive regime (see Figures and b). The obtained 
values for D are compared to the ones expected from the 
diffusion approximation assuming Dda = vlt/3. Figure 
shows a subset of simulated mean square width slopes 
obtained from samples with It = L. 

As can be seen, high deviations are obviously to be 
found at such a low optical thickness, especially in the 
proximity of n = 1, which is unexpected since this con¬ 
dition is sometimes quoted as a safer configuration for 
the diffusion approximation M- It is also worth noting 
that the mean square width slope of simulated data is 
always, over the whole set of parameters [12], equal to or 
greater than the value expected by standard diffusion with 
the ratio approaching 1 for increasing optical thicknesses, 
thus recovering the diffusive prediction as expected. This 
general trend can be intuitively explained considering the 
d-dimensional modeling of diffusion as a random walk pro¬ 
cess, which, given any arbitrary step length distribution 
P{1) with finite moments (1) and (P), predicts a mean 
square d-dimensional displacement growing as 2dDt with 


D = 


1 


2d (0 



( 2 ) 


where the last equality holds for an exponential step length 
distribution with average step length It |S] . Indeed, as the 
optical thickness of the simulated slab decreases, trans¬ 
port occurs in an increasingly planar geometry. Hence, 
as suggested by eq. ([^, the effective diffusion coefficient 
D as inferred from the mean square width slope might 
be up to 3/2 times higher than its bulk nominal value. 
The perceived environment dimensionality is also clearly 
affected by the refractive index contrast. Near n = 1, 
where deviations are more relevant, any late-time surviv¬ 
ing energy packet will have performed an almost planar 
trajectory, more similar to a purely 2-dimensional walk. 
On the contrary, with increased reflections at boundaries. 
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Figure 1. Sketch of the investigated configuration. An infinite slab is illuminated by a pencil beam 5{t) pulse. Transmitted light 
is collected at different times and positions. A few (longer) trajectories and normalized transmitted intensities are presented 
for illustrative purposes in the case of an optically thin, index matched slab with thickness and scattering mean free path 
equal to 1mm (i.e. optical thickness = 1), shown to scale. The scattering anisotropy g is set to 0. An approximately Gaussian 
profile is transmitted at each time slice, whose mean square width and integrated intensity respectively grows linearly and falls 
exponentially. 


trajectories are allowed to fold back into the sample thus 
perceiving a more 3-dimensional environment. 

A closer look on the data reveals an unexpected feature 
occurring in the close proximity of n = 1. A set of g- 
cross-cuts in the plane It/L^ = 1 (Figure]^) reveals a 
sharp modulation of the mean square width slope across 
unitary index contrast. In particular, contrarily to what 
one would expect from the above discussion, the mean 
square width expansion exhibits a local minimum at n = 1 
rather than a maximum. In its close proximity, diffusion is 
asymmetrically enhanced reaching an absolute maximum 
around n = 1.016 for g = 0. These features arise from 
a subtle interplay between the optical properties of the 
sample, as we will discuss in the following section. 

Before any further discussion on the mean square width 
expansion, it is anyway worth noting that a qualitatively 
similar behavior is mirrored in the lifetime relative de¬ 
viations from the diffusive prediction, and as such this 
effect is not strictly limited to transverse transport. This 
is of particular interest since lifetime measurements have 
long been experimentally accessible and exploited to esti¬ 
mate the diffusion coefficient via the relation (in case of 
a non-absorbing medium) 


where Leg = L + 2ze is the effective thickness of the 
medium and Ze represents the extrapolated boundary 
conditions for a given refractive index contrast. The simi¬ 
larity in the dependence on n is particularly evident when 
plotting the same cross-cuts as for the mean square width 
slope relative deviations (Figure [^). Two main features 
are worth commenting. Firstly, it is evident that the life¬ 
time ratio can deviate more significantly from the diffusion 
approximation than the mean square width slope. Sub¬ 
stantial deviations from the lifetime prediction are indeed 


known to persist also at higher optical thicknesses, espe¬ 
cially for n > 1 despite the correct extrapolated bound¬ 
ary conditions being applied II1I5D]. As shown by our 
complete set of simulations m, the mean square width 
prediction error at OT = 10 does not exceed « 2 % for 
each refractive index contrast and scattering anisotropy 
probed, while the same deviations for the lifetime can 
exceed 20% for n > 1. The second relevant difference 
regards the fact that the t/tda ratio can evidently take 
values both above and below 1, depending subtly on the 
scattering anisotropy and the refractive index contrast 
of the sample. This observation might provide a clearer 
picture on why a diffusion-based retrieval of the diffusion 
coefficient from a lifetime measurement is sometimes re¬ 
garded as a poor estimation, since it can lead both to over- 
or underestimated values (Figure]^). This is further il¬ 
lustrated in Figure for a couple of representative cases 
exhibiting opposite deviations that can persist even at 
higher optical thicknesses. To date indeed, experimental 
data and theoretical predictions are inconsistent. While 
the former bring generally evidences suggesting a decreas¬ 
ing diffusion coefficient with decreasing thickness mull], 
the latter have so far mainly provided arguments in favor 
of the opposite behavior [niiniissiiii]- In this respect 
it is worth noting that our simulations suggest that there 
is a region in the parameter space where the t/tda ratio 
could exceed 1, which can lead to the observed decreas¬ 
ing diffusion coefficient with decreasing thickness. The 
analysis on the decay lifetimes confirms the importance 
of an accurate and precise modeling of the index contrast, 
which we think has been often overlooked, for example 
by considering a symmetric averaged contrast to model 
asymmetric experimental configurations [n [201 Ha¬ 
lt is worth noting that a better modeling of the T{lt) 
dependence can be envisioned, given its inappropriate¬ 
ness at low optical thicknesses. In the standard diffusion 
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Figure 2. (a) Cross-cuts along different values of g of the 
simulated hyper-surface dataset m, taken a,t It = Lq = 1 mm 
(OT = 1). For each point, the MSW slope has been evaluated 
from a simulation of 10^° energy packets. The introduction 
of a small amount of boundary reflection appears to enhance 
in-plane diffusion remarkably, which also affects g > 0 and 
OT > 1 samples m- A subset of parameters (black filled 
dots) where the effect is more dramatic is further investigated 
in the following, based on a much larger statistics, (b) MSW 
growth for the three different configurations highlighted in the 
upper panel. Linear fits are performed excluding t < 4r. 


approximation, the lifetime equation can be written as 


{L + 2Ze)^ 
tt'^D 


— Po + P+ih + P-ih ^ 


(4) 


where po = SALjn'^v, = 16A^/37r^n, p_i = 
and A = A{n) is the correction factor accounting for 
internal reflections. By including also higher order terms 
as p+2^t ^'^'4 it is indeed possible to perfectly 

reproduce each simulated value (Fig. [^, red curves) on 
a broad parameter space, while preserving important 
physical properties of the lifetime dependence such as its 
divergence with OT —> oo, which would not be guaranteed 
by a generic polynomial fit. 


Figure 3. (a) Detailed view of the n dependence of the lifetime 
deviation at OT = 1. Especially for lower g factors, the ratio 
can clearly go above 1 for specific contrasts. The dependence of 
T on OT for values highlighted as black filled circles is shown in 
panel (b), along with the DA prediction and a higher order fit. 


IV. DISCUSSION 


In order to explain the origin of the features shown 
above for the mean square width dependence on n, we 
now focus on three significant sets of parameters (high¬ 
lighted as black filled circles in Figure]^) representing 
key points of the observed peak, i.e. n = 1,1.016 and 1.1, 
the last being the coordinate where diffusion recovers the 
value in n = 1. These three particular configurations were 
further investigated separately with the aim of collecting 
detailed statistics at long times, with lO^'^, 0.5 x 10^^ and 
10^^ energy packets each. The unprecedented magnitude 
of these simulations and the subtlety of the effects investi¬ 
gated required the development of dedicated software with 
a particular focus on numerical stability, precision and 
reproducibility, as described elsewhere m- In particular, 
dealing with more than « 10^*^ exponentially distributed 
random variates generally requires the use of Pseudo Ran¬ 
dom Number Generators (PRNGs) with more than 32 
random bits and sufficient floating point precision (i.e. 
long double) in order not to introduce statistically signifi¬ 
cant truncations in the drawn values, both of which are 
features not readily available in existing MC solutions. 
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As suggested by eq. ([^, the most straightforward in¬ 
sight on the effective D is obtained by directly looking 
at the distribution of the step lengths performed during 
the random walk. In principle, each energy packet is 
simulated according to the same step length distribution 
P{1) = exp—l/ls, which is characterized by a decaying 
slope and an average value of 4- On the contrary, we 
found that there exists a clear correlation, induced by 
this confined geometry, between a long permanence inside 
the sample and an unevenly sampled step length distribu¬ 
tion. Figure]^ shows the histograms of the step lengths 
and scattering angles between two consecutive scattering 
events of energy packets transmitted at t = t/run = 90 ps 
(corresponding to a path length of « 27Lo) compared to 
their nominal distributions (dashed blue lines). 

The step length distributions exhibit enhanced tails in 
this confined geometry for all the three simulated refrac¬ 
tive index contrasts. This is consistent with a diffusion co¬ 
efficient higher than the one expected by the diffusion ap¬ 
proximation, because heavier tails will generally be associ¬ 
ated with a greater enhancement of the second rather than 
the first moment (cfr. eq. ([^). The peculiar late trans¬ 
port regime arising in this confined geometry causes the 
nominal step length distribution provided by the PRNG 
to be sampled unevenly: despite the fact that a long step 
in a very thin sample will generally cause the packet to 
exit the slab, those few that happen to remain inside will 
be able to reach long surviving times without undergoing 
many scattering events. Notably, in the case of refractive 
index contrasts very close to 1, the distribution of the 
step lengths features a selective enhancement of the longer 
values, which is slightly more marked for n = 1.016. This 
might be due to the fact that, for such a small refractive in¬ 
dex contrast, total internal reflection is already significant 
{9c = 79.8°). Consider as an example the case of an energy 
packet taking a very long step: if internal reflections are 
absent, extremely narrow angular conditions must hold in 
order for it not to exit the slab. Conversely, even a tiny 
contrast allows to largely relax such condition, introduc¬ 
ing a significant increase in the survival probability of a 
long-stepping energy packet while only marginally affect¬ 
ing others. In other words there’s a positive correlation 
between long steps and shallow incidence angles, whose ef¬ 
fects become more apparent when such angles are the only 
ones undergoing total internal reflection (which also ex¬ 
plains why the enhancement is asymmetric around n = 1). 
On the other hand, with increasing contrast, more energy 
packets will be held inside the slab irrespective of their 
incidence angle (and hence of the length of their step), 
thus washing out the observed MSW spreading boost. 

Interestingly, the sampling of the angular variables is 
also modified at late times, as shown in Figures |^-c for 
the same set of simulations. While tracing each energy 
packet trajectory, the cosine of the scattering (polar) an¬ 
gle 9 and the azimuthal angle ip are generated through the 
PRNC so to be uniformly distributed in [—1,1] and [0, 27r), 
respectively. On the contrary, the cos 9 distribution at late 
times exhibits a back and a forward peak, with a plateau 


around cos 9 = 0. Further insight on these features is 
gained by looking at the concurrent modifications in the 
Ip statistics, which exhibit two symmetric peaks around 
■P-K12. These two peaks correspond to a right/left turn in 
the slab plane, which helps keeping the trajectory inside 
the sample irrespective of the scattering angle 9. The two 
cos 9 peaks can also be intuitively understood by consid¬ 
ering a typical step in a very long trajectory: this will 
generally be a long step (i.e. / ^ 4 = Lq) mostly aligned 
with the slab plane. As such, scattering angles close to 
0 = 0° or 180° will guarantee that the trajectory will 
continue within the slab irrespective of what azimuthal 
angle is drawn. Actually, since a typical step will not be 
in general perfectly parallel to the interfaces, a scattering 
angle of 0 « 180° should provide higher chances of staying 
inside the sample, hence its higher probability. Interest¬ 
ingly, this results in a cos 9 distribution with a slightly 
negative average value, which is also able to influence the 
effective diffusion properties exhibited by the sample. 

The time dependence of both {cos 9) and the ratio 
(P)/2{l) is tracked in Figure]^ along with their nomi¬ 
nal values (dashed blue lines). Each point collects the 
statistics of those energy packets that were transmitted 
within that time bin. The obtained curves successfully 
validate our interpretation based on a random walk pic¬ 
ture of the diffusive process as expressed by equation 
([^, exhibiting good qualitative agreement with the mean 
square width expansion data. In principle, the overall 
observed diffusion process will be influenced by both the 
modified step length and angular statistics, which in the 
investigated configurations appear to have opposite effects. 
While the latter would indeed tend to slightly slow down 
diffusion, the predominant effect is clearly coming from 
the step lengths being substantially increased, leading to 
the observed enhanced in-plane diffusion especially for 
n = 1.016. Notably, different configurations might lead 
to a different overall balance between these two effects, 
which also appear to saturate to their respective asymp¬ 
totic values with slightly different time scales, further 
illustrating the need for additional investigation even for 
the simple homogeneous and isotropic single slab model. 

These results show that the currently established the¬ 
oretical framework linking radiative transfer theory to 
diffusion needs to be further refined, especially for thin 
systems and more generally for confined geometries. In 
particular, concerning microscopic optical properties such 
as the scattering anisotropy or the mean free path, it 
seems appropriate to introduce a distinction between an 
intrinsic and an effective counterpart, where the former 
is the one that we are typically interested in retrieving 
while the latter might have a very different value and 
nature (e.g. tensorial instead of scalar) depending on in¬ 
cidental geometric conditions. We have shown in fact 
that a homogeneous, isotropic slice of a certain disordered 
material will eventually reach an asymptotic, multiple 
scattering regime characterized by different (and possibly 
anisotropic) statistical distributions, depending on both 
sample size and boundary conditions. 
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(a) 



Figure 4. Late-time modification of the step length and angular distributions. Three simulations with OT = 1, y = 0 and 
n = 1,1.016 and 1.1 are considered. Panel (a) shows the histogram of the distribution of step lengths between consecutive 
scattering events performed by energy packets transmitted at t = 90 ps. Also angles become unevenly sampled at late times, as 
shown in panels (b) and (c). Each blue dashed line represent the actual distribution fed into the simulation. 


It should be noted that such discrepancy — which we 
are now able to correctly identify as the emerging of an 
effective transverse diffusion coefficient — has already 
been reported experimentally in samples with an optical 
thickness as high as 8 [TH], suggesting that it can indeed 
represent an appreciable issue in a broad range of appli¬ 
cations. It is worth stressing again the different nature of 
mean square width and lifetime deviations at higher opti¬ 
cal thicknesses. Even at OT = 25, for a slab sample with 
a typical glassy refractive index contrast, the diffusive 
prediction for the lifetime is still wrong by more than 2 % 
while the mean square width expansion discrepancy is one 
order of magnitude smaller, suggesting that the lifetime 
deviation is due to a defective diffusion-based modeling of 
boundary conditions rather than to an actually different 
effective diffusion coefficient. 

The asymptotic nature of the effective diffusion coeffi¬ 
cient in a thin slab is further illustrated in Figure]^ where 
the time evolution of the step length distribution is shown 
for n = 1 (the n = 1.016 and 1.1 cases are analogous). 
The time-resolved distributions seem to converge towards 
a single asymptotic envelope distribution with a well 
defined asymptotic decay which seems to be uniquely de¬ 
termined by the properties of the sample. It is interesting 
to compare the histogram of the actual steps performed 
inside the sample with the histogram of the ones drawn 
from the PRNG. The two differ only for the last step, 
whose length is respectively considered either partially 
(up to the interface) or totally. As can be seen, only the 
earliest transmitted energy packets did actually sample 
the expected distribution evenly: in this respect, while 
on one hand they underwent too few scattering events 
and thus fall outside the diffusive approximation, on the 
other hand they are the only ones to migrate according 
to the nominal step length distribution. 

Even if at an optical thickness close to 1 it is indeed 
questionable to expect diffusion predictions to hold, it is 
nonetheless true that we were able to follow the evolu¬ 


tion of light deeply into the multiple scattering regime, 
where light is undoubtedly diffusing as the radial inten¬ 
sity profiles and their linear variance expansion clearly 
demonstrate. Following the usual approach based on the 
diffusion approximation, it may be argued that it could 
be possible to take into account all the effects discussed 
here by means of some refined extrapolated boundary 
condition. On the contrary, the optical thickness must be 
regarded as one of the parameters actively and indepen¬ 
dently affecting transport properties. Moreover, extrapo¬ 
lated boundary conditions are devised to impact and sub¬ 
stantially correct quantities such as total transmittance 
which, conversely, would be hardly affected by asymptotic 
modifications of the effective diffusion coefficient. 

The absolutely non-trivial dependence of the random 
walk statistics near n = 1 should also be stressed. This 
is particularly relevant since index matched boundary 
conditions are customarily considered as a “forgiving case” 
where the diffusive approximation is expected to hold 
better also for thinner systems [14]. While this might be 
true under various circumstances, from an experimental 
point of view it is clear that typical attempts to approx¬ 
imately index-match a thin sample within a reference 
material might more likely introduce an error rather than 
neutralize it, since it’s the very smallest index mismatch 
that results in the strongest deviations from the theory. 


V. CONCLUSIONS 

We have studied the diffusive transport regime in opti¬ 
cally thin samples, which develops at late times mainly 
within the slab plane and shows peculiar features espe¬ 
cially for small refractive index contrasts. In order to 
explain our numerical results we performed extensive sim¬ 
ulations of three key configurations which allowed us to 
inspect transport properties on exceptionally long time 
scales, so far unexplored. Our investigation revealed a 
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Figure 5. Time evolution of the effective scattering length 
and effective scattering coefficient as inferred from the time- 
resolved effective step length and angular distributions. For 
each time bin, (cosS) and {l^)/2{l) are calculated from the 
retrieved time-dependent distributions, based on the trajec¬ 
tories of energy packets transmitted within that time bin. 
Blue dashed lines represent the expected values for the two 
distributions. 



Figure 6. Time evolution of the step length distribu¬ 
tion for n — 1 for energy packets transmitted at f = 
10,20,30,40,50,60,70,80 and 90 ps. Black and gray curves 
show respectively the histograms of the lengths of the steps 
taken inside the sample and drawn through the PRNG. The 
two differ for the last step, whose length is accounted for 
only up to the exit surface in the first histogram. At late 
times the two sets of curves become indistinguishable since, 
as expected, the contribution of the last step to the whole 
trajectory becomes eventually statistically negligible. 

subtle interplay occurring between the actual thickness 
of the slab, the refractive index contrast and the scatter¬ 
ing anisotropy. This results in a regime that is properly 
diffusive but which cannot be described in terms of the 
standard diffusive approximation. Such interplay gives 
rise to a different diffusion coefficient which emerges nat¬ 
urally from the overall optical and geometric boundary 
conditions of the sample, and is univocally determined 
by them through yet unknown relations. In this respect, 
our findings recall a recently published work where it is 
analogously demonstrated that the link between micro¬ 


scopic (i.e. the scattering coefficient) and macroscopic (i.e. 
the diffusion coefficient) transport parameters remains 
unknown for diffusive anisotropic media [S], therefore re¬ 
inforcing the claim for a rethink of our present picture of 
the diffusive modeling of a random walk-based transport 
problem. 

As we have shown in Figures and transport eventu¬ 
ally converges to an asymptotic regime characterized by a 
well-dehned diffusion coefficient. Interestingly, we found 
that because the migration of simulated energy packets 
is directly altered by the sample configuration, the sam¬ 
ple itself appears generally to be less scattering than it 
actually is. In other words, once the diffusive regime 
is reached, energy packets migrate as if scatterers were 
further apart than they actually are, i.e. with an effective 
transport mean free path which is greater than the one 
intrinsic to the material. From this perspective, we have 
demonstrated that even in the presence of clear diffusion 
signatures (such as a steadily linear mean square width 
growth, see Fig. [^) care has to be taken on whether it is 
appropriate to resort to standard diffusive approximation. 

Because of the asymptotic nature of the effects here 
described, when studying thin samples only a small frac¬ 
tion of incoming light is actually subject to this effective 
transport mean free path. Nonetheless experimental tech¬ 
niques capable of detecting this discrepancy have been 
recently reported |18j . Furthermore other applications 
can be envisioned where multiple scattering in thin lay¬ 
ers, even if limited to few energy packets, could play a 
significant role (e.g. random lasers). 

Finally, despite our case study is in the field of light 
transport, random walks are extremely general models for 
a broad range of phenomena in complex systems, from 
molecular kinetics to social behaviors. Therefore our 
findings could have far-reaching consequences in other 
scenarios where tight geometric confinement holds and 
analogous boundary conditions can be modeled. In this 
respect we demonstrated how, depending on the applica¬ 
tion, the interplay between transport properties and the 
environment geometry can give rise to sharp and unex¬ 
pected macroscopic migration features, which is either a 
condition to be taken into account or one to be possibly 
exploited as an engineering degree of freedom to enhance 
transport. 
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